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Simultaneous regulation of cell size and chromosome replication in bacteria 
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Bacteria are able to maintain a narrow distribution of cell sizes by regulating the timing of cell 
divisions. In rich nutrient conditions, cells divide much faster than their chromosomes replicate. This 
implies that cells maintain multiple rounds of chromosome replication per cell division by regulating 
the timing of chromosome replications. Here, we show that both cell size and chromosome replication 
may be simultaneously regulated by the long-standing initiator accumulation strategy. The strategy 
proposes that initiators are produced in proportion to the volume increase and is accumulated at 
each origin of replication, and chromosome replication is initiated when a critical amount per origin 
has accumulated. We show that this model maps to the incremental model of size control, which 
was previously shown to reproduce experimentally observed correlations between various events in 
the cell cycle and explains the exponential dependence of cell size on the growth rate of the cell. 
Furthermore, we show that this model also leads to the efficient regulation of the timing of initiation 
and the number of origins consistent with existing experimental results. 


INTRODUCTION 

Bacterial cells are extremely proficient in regulating 
and coordinating the different processes of the cell cy¬ 
cle. The Cooper-Helmstetter model proposes a molecular 
mechanism that couples two such processes, the replica¬ 
tion of the chromosome and the division of the cell |T] . In 
the model, cell division occurs a constant duration after 
the initiation of chromosome replication. The model im¬ 
plies a tight coordination between replication initiation 
and cell division such that in cells able to double faster 
than their chromosomes can replicate, multiple rounds of 
replications proceed simultaneously m- To answer how 
cells regulate the timing of initiation, it was proposed 
that replication initiation factors accumulate to a criti¬ 
cal amount per origin of replication to trigger the initia¬ 
tion of replication [3]. Since the conception of the above 
model, many experiments and models have attempted 
to capture the molecular mechanisms responsible for the 
initiation of multiple rounds of replication. However, no 
model has been completely satisfactory j3j. 

As a result of the coupling between replication and 
division, the average cell size per origin is approximately 
a constant independent of the growth rate of the cell [5]. 
Furthermore, it is now understood that a common size 
regulation strategy for organisms including bacteria and 
budding yeast is the incremental model in which division 
occurs upon the addition of a constant size dependent on 
the growth rate of the cell EHS|. However, the molecular 
mechanisms responsible for the incremental model of size 
control remain in question. 

Our main result in this work is to show that the ini¬ 
tiator accumulation strategy not only regulates size ac¬ 
cording to the incremental model, but also regulates si¬ 
multaneously the timing of initiation and the number of 
origins of replication. The strategy says that replication 
initiates upon the accumulation of a critical amount of 
replication initiation factors per origin. We emphasize 
the importance of the partitioning of replication initia¬ 


tion factors amongst origins, which we show is essential in 
order for the multiple replication forks to be adequately 
regulated. We show, analytically and numerically, that 
this strategy robustly regulates both cell size and the 
number of origins. Agreement between existing exper¬ 
iments and predictions of the above model reveals es¬ 
sential features that must be captured in any molecular 
mechanisms coordinating replication initiation and cell 
division. Finally, we make distinct predictions regarding 
the distribution of cell sizes at initiation of replication. 


METHODS 

Multiple origins accumulation model 

We consider the regulation strategy in which replica¬ 
tion initiates upon the accumulation of a critical amount 
of replication initiation factors, or ’’initiators”, per ori¬ 
gin of replication [3j. We assume that the initiators are 
expressed via an autorepressor model, as seen in Fig¬ 
ure mum. In this model, a protein is expressed such 
that its concentration c remains constant and indepen¬ 
dent of the growth rate of the cell, which is plausible to 
achieve through autorepression. Therefore, an increase 
in the volume of the cell corresponds to a proportional 
increase in the copy number of this autorepressing pro¬ 
tein. A second protein is the initiator and is expressed 
under the same promoter as the first, but in contrast to 
the first protein, it is localized at the origins of replica¬ 
tion. For simplicity, we assume that the initiators are 
equally partitioned amongst the origins. Initiation then 
occurs when a critical copy number per origin Writicai of 
the localized initiators is reached, after which the initia¬ 
tors are assumed to degrade. Under these assumptions, 
the copy number of the initiator effectively measures the 
increase in volume since initiation. 

More precisely, if a cell initiated a round of replication 
at volume V{ into O number of origins, the amount of ini- 
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tiators initiators immediately after initiation is zero. To 
initiate the next round of replication, the cell must accu¬ 
mulate O-iriticai initiators, but because the initiator is 
expressed under the same promoter as the autorepressor, 
the cell must also accumulate Oiriticai autorepressors. 
Because the concentration of the autorepressor is con¬ 
stant, this implies that the cell must accumulate a corre¬ 
sponding volume A = -/V cr iticai/c per origin, independent 
of the growth rate, to trigger the next initiation. 

Thus, on a phenomenological level, the above biophys¬ 
ical model maps to the following regulation strategy for 
initiation, 


tot,next i /n a /i\ 

Vi &Vi + O A, (1) 

Eq. Q says that if a cell initiated a round of replication 
at cell volume Vi into O number of origins, then the cell 
will attempt to initiate another round of replication at 
total volume ^ ot > next ) which is the sum of the volumes 
of all cells in the lineage since the initiation event at u, 
(typically two cells). This is not to be confused with the 
threshold model in which cells initiate upon reaching a 
threshold volume proportional to the number of origins, 
uf ext oc O. For the rest of this manuscript, O will denote 
the number of origins after initiation at cell volume Vi 
but before initiation at total cell volume v * ot ’ nex \ 

We assume an exponential mode of growth for cell vol¬ 
ume with a constant doubling time r and a corresponding 
constant growth rate A = ln2/r [11] . From Eq. [I] and 
the exponential mode of growth, durations between ini¬ 
tiations are 


1, ( O A\ „ 

U ~ A hl ( 1 + Vj + ^’ 


( 2 ) 


where f represents some noise in the initiation process. 
An initiation event will trigger a division event after a 
constant duration C+D , where C and D are respectively 
the constant duration required to replicate the chromo¬ 
some and the constant duration between replication ter¬ 
mination and division pTj. We will refer to Eq. |T| as the 
multiple origins accumulation model (i.e. initiators are 
accumulated per origin). Figure [2] illustrates this regula¬ 
tion strategy. We note that the strategy described here 
is mathematically equivalent to the ’’replisome” model of 
Bleecken jT2] (not to be confused with the current use of 
the term replisome). 

Finally, we will not take into account additional bio¬ 
logical mechanisms that act at the level of the initiation 
of chromosome replication, such as oriC sequestration, 
Dam metliylation, and the ’’eclipse” phenomenon d- 
115] . While these mechanisms are important to prevent 
rapid re-initiations, by themselves they are insufficient in 
ensuring an appropriately coordinated coupling between 
chromosome replication and cell division, which is the 
main focus of our work. 
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FIG. 1. Schematic of the autorepressor model of initiator 
expression. An autorepressor is expressed such that its con¬ 
centration remains constant and independent of the growth 
rate of the cell. The initiator is expressed in proportion to 
the autorepressor, but is localized at the origins of replica¬ 
tion. Adapted from Sompayrac and Maaloe [El- 
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FIG. 2. Schematic of the regulation strategy of the multiple 
origins accumulation model. See text for the details of the 
model. Slow growth denotes 0 < C+D < 1. Faster growth 
denotes 1 < °+ D . In the above example, C+D < 2. 


Numerical simulations 

We can numerically simulate the multiple origins accu¬ 
mulation model given C + D, r, and A as experimentally 
measurable parameters. First, we initialize a population 
of N cells with uniformly distributed cell ages. Durations 
between initiations are calculated as Eq. [2] and the noise 
in the initiation process is assumed to be normally dis¬ 
tributed with standard deviation a Tl though the precise 
nature of the noise does not affect any of our conclusions. 
It is assumed that in an initiation event, the number of 
origins in a cell is doubled. The corresponding division 
event occurs after a constant time C + D. In a division 
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Generations 


FIG. 3. ti, O, and Vb approaching stationary distributions 
in numerical simulations of the multiple origins accumulation 
model. First, we initialize a population of N cells with uni¬ 
formly distributed cell ages. Durations between initiations 
of replication are calculated as Eq. [2] and the noise in the 
initiation process is assumed to be normally distributed. In 
an initiation event, the number of origins in a cell is doubled. 
The corresponding division event occurs after a constant time 
C + D. In a division event, the number of origins in a cell, 
along with the size of the cell, is halved, and two identical 
cells are generated. There are no division events without the 
corresponding initiation events. Following this procedure, a 
population of cells will robustly reach a stationary distribu¬ 
tion of cell sizes and number of origins per cell regardless of 
initial conditions. The plots here track one lineage of cells. 
Here, C+D = 70 mins, r = 20 mins, and ov = 4 mins. These 
are biologically realistic choices. We set A = l/2 ( - c ' +i W’ r so 
that ( Vb) ~ 1. 


event, the number of origins in a cell, along with the size 
of the cell, is halved, and two identical cells are gener¬ 
ated. We neglect the stochasticity arising from asymmet¬ 
ric divisions, which do not significantly affect any of the 
results. There are no division events without the cor¬ 
responding initiation events. Following this procedure, 
a population of cells will robustly reach stationarity re¬ 
gardless of initial conditions, as seen in Figures [3] and 

m 


RESULTS 

Multiple origins accumulation robustly and 
efficiently regulates the number of origins of 
replication 

An important measurable consequence of the tight cou¬ 
pling between replication initiation and cell division is 
the average number of origins of replication per cell. It 
has been theoretically shown that the average number of 
origins per cell is 

(O) = 2^ C+D ^ T . (3) 



FIG. 4. Stationary exponential distribution of cell ages. 
Simulations are the same as Figure [3] The line plots p (a) = 
(ln2)2 1_ “, as described in the Appendix. 

The derivations leading up to Eq. [3j summarized in Ap¬ 
pendix, hinge on assuming an efficient process: the pop¬ 
ulation is growing exponentially and reaches a stationary 
distribution of cell ages, implying that there are no de¬ 
lays due to chromosome replication. These assumptions 
are independent of any molecular mechanisms for initia¬ 
tion and thus should be fulfilled by any efficient mecha¬ 
nism. We will show that the multiple origins accumula¬ 
tion model reproduces Eq. [3] 

The multiple origins accumulation model regulates ini¬ 
tiation via negative feedback on the volume at initiation 
described by Eqs. [l] and [2] The feedback enables cells to 
maintain a stationary average volume at initiation and 
a stationary average duration between initiations despite 
noise in the initiation process. Specifically, if a cell ini¬ 
tiated replication at volume Vi = OA, then the duration 
to the next initiation event is 

< 4 > 

But if a cell initiated replication at a slightly larger vol¬ 
ume Vi = O A + 5v, the duration to the next initiation 
event is 

^ 0 > 4‘"( 1 + ra ) ST ' < 5 > 

Eqs. [4] and [5] say that a cell that initiated at a slightly 
larger volume than average tend to initiate again faster 
than average so that its volume at next initiation is again 
near the average. Similar reasoning says that cells that 
initiated at slightly smaller volumes tend to initiate again 
slower than average. In this way, cells maintain a station¬ 
ary average volume at initiation and a stationary average 
duration between initiations. 
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Furthermore, the feedback enables cells to maintain a 
balanced cell cycle, in which there is on average one and 
only one initiation event per cell cycle. In the case of 
negligible noise, a balanced cell cycle implies that cells 
will initiate at cell age uni 


= 1 + L 


G 


Dj C 


D 


( 6 ) 


where a = 0 represents cell birth, a = 1 represents 
cell division, and |J is the mathematical floor operator 
(largest integer smaller or equal to the argument). But 
in the case of realistic noise, a cell may initiate an ex¬ 
tra round of replication if the noise is negative enough, 
£/r 1$ [{C+D)/t\ — (C + D)/t , which corresponds to an 
extra initiation at volume v[ = 20 A — Sv. The multiple 
origins accumulation model is robust to these stochas¬ 
tic events because a cell that initiated an extra round of 
replication will initiate again after 



< 7) 

In other words, cells with extra rounds of replication will 
initiate slower than those without so that the station¬ 
ary average duration between initiations is maintained. 
The cell cycle following the extra initiation will typically 
not have any initiations, so that the initiation following 
the extra initiation will occur at approximately the aver¬ 
age volume at initiation. A cell that missed a round of 
replication will return to a balanced initiation process in 
the analogous manner. In this way, the multiple origins 
accumulation model is able to efficiently maintain a bal¬ 
anced cell cycle in fast growth conditions. In contrast, 
the model simulated by Campos et al. [9| is not robust to 
the noise in the initiation process, because in their model, 
the incremental volume needed to trigger initiation is not 
partitioned between origins. We will elaborate on this in 
Section 3.3. 

The multiple origins accumulation model is therefore 
able to robustly regulate the timing of initiations in face 
of extra initiations. Extra initiations can occur not only 
because of noise in the initiation process, but also be¬ 
cause of a shift in the growth rate of the cell such as that 
found in a shift-up experiment, in which a population of 
cells is abruptly switched from one nutrient condition to a 
richer nutrient condition allowing for faster growth. The 
increase in growth rate corresponds to a decrease in the 
duration between initiations. Therefore, cell in a shift-up 
experiment will initiate extra rounds of replication in the 
cycle immediately following the shift-up, but as we have 
seen, the multiple origins accumulation model is able to 
appropriately regulate the timing of initiations to reflect 
the new growth rate. Simulations of the multiple origins 
accumulation model reached stationary distributions of 
cell ages, durations between initiations, cell sizes, and 
number of origins per cell, regardless of initial conditions 


FIG. 5. ( O) as a function of r shows regulation of chromo¬ 

some replication. Simulations are the same as Figure [3] with 
a varying r = 20 to 100 mins, a fixed C + D = 70 mins, and 
oy/r = 0.2. Dashed line plots Eq. [3] Similarly, inset plots 
log (O) as a function of A 


or the magnitude of the noise £ in Eq. [2j as seen in Fig¬ 
ures |3] and |4] Simulations also show that the number of 
origins is regulated as in Eq. [3j as seen in Figure [5] The 
above considerations show that the multiple origins accu¬ 
mulation model regulates the number of origins robustly 
and efficiently in face of noise in the initiation process. 

Multiple origins accumulation robustly regulates cell 

size 

It was recently shown that the multiple origins accu¬ 
mulation model of replication initiation reduces to the 
incremental model of size regulation [B] 

v d ttv b + v 0 . ( 8 ) 

Eq. [8] says that if a cell is born with volume v b , then 
the cell will attempt to divide at volume v d , where vq is 
the constant incremental volume from birth to division. 
In fact, Vo can be expressed in terms of known param¬ 
eters. First, if a cell initiated replication at volume i>i, 
then the corresponding division event will occur at to¬ 
tal volume Vj ot = Vi2^ c+D ^ T . But there will have been 
log 2 0 division events since initiation at Vi, so that the 
corresponding volume at birth is 

_ v^ ot _ v. 1 2 {C+D)It 


In a balanced cell cycle, the next initiation event will 
occur at total volume ^ot.next ~ Vi -g QA and the 
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corresponding division event will occur at total volume 
^tot,next _ ^tot,next 2 (c+D)/r^ Similarly, there will have 

been log 2 O division events since initiation at u" ext , so 
that the corresponding volume at division is 


Vd 


v 


tot, next 
d 


o 


Vi 2 {C+D)/T 

O 


+ A2 i C + D )/T. 


( 10 ) 


Therefore 


v d ~v b ^ A2^ C+D ^ T . ( 11 ) 


Within the multiple origins accumulation model, this 
derivation is valid for any C + D and r [ 6 ] . 

The incremental model of size regulation predicts dis¬ 
tributions, correlations, correlation coefficients, and scal¬ 
ings consistent with existing measurements PH- In par¬ 
ticular, the average cell volume at birth 

(v b ) » A2^ c+d ^ t . ( 12 ) 


Eq. [l2]says that the average cell volume at birth is expo¬ 
nentially dependent on the growth rate, a well-known and 
well-tested result for E. coli and B. subtilis 8 , 17, 18 . 
Simulations of the multiple origins accumulation model 
also confirm this result, as seen in Figure [7} Thus, the 
multiple origins accumulation model robustly regulates 
cell size. 


Master accumulation predictions are inconsistent 
with existing experiments 


Consider the regulation strategy 

tot,next _ , a 

Vi ~ Vi + A. 


(13) 


Eq. 13 says that if a cell initiated a round of chromosome 
replication at cell volume Vi with O number of origins of 
replication, then the cell will attempt to initiate another 
round of replication at total volume z ,* ot ’ next ; w here A is a 
constant volume independent of the growth rate. This is 
an incremental model of size control applied at initiation. 
If we assume the same mode of initiator expression as 
before, then the regulation strategy described by Eq. [13] 
corresponds to replication initiation upon accumulation 
of a critical amount of initiators without partitioning of 
initiators between origins. Instead, a plausible molecular 
picture is that of initiators accumulating at a ’’master” 
origin, whose initiation triggers the cascade initiation of 
other origins [T5j. We will therefore refer to Eq. 13 as 
the master accumulation model. In contrast to the mul¬ 
tiple origins accumulation model described above, here 
in the master accumulation model, the total volume at 
next initiation does not depend on the number of origins 
present in the cell. 

As before, we assume an exponential mode of growth 
with a constant doubling time r. The durations between 


initiations are therefore 


t% (O) — In ( 1 + — ) + Ci 


(14) 


where £ represents some noise in the initiation process. 
Again as before, an initiation event will trigger a division 
event after a constant duration C+D. Eq. [14| differs from 
Eq. [2] by a missing factor of O. As we show below, the 
factor of O is essential in regulating appropriately the 
timing of initiations and the master accumulation model 
does not reproduce the well-known exponential scaling 
of cell size with growth rate. The derivation follows and 
is similar to that of the multiple origins accumulation 
model. First, if a cell initiated replication at volume 
Vi, then the corresponding division event will occur at 
total volume v ^ ot = Vi2 ( ' C+D ' ) / T . But there will have been 
log 2 0 division events since initiation at Vi, so that the 
corresponding volume at birth is 


v^ _ v i 2( c+D '>/ T 

~o ~ O 


(15) 


In a balanced cell cycle, the next initiation event will 
occur at total volume ^ ot > next ~ Vi 4 . A and the cor¬ 
responding division event will occur at total volume 
^tot,next _ ^tot,next 2 (c+D)/r_ Similarly, there will have 

been log 2 O division events since initiation at i>" ext , so 
that the corresponding volume at division is 


v d 


v 


tot, next 
d 


o 


v>2 {C+D)/t 

O 


A 2 (C+D)/t 

O 


(16) 


Therefore, 


A2 (C+D)/r 

v d -v b = ---. (17) 

This derivation is valid for any C + D and r. However, 
from Eq. [3j O should scale exponentially with the growth 
rate like 2S C+D )/ T so that 


( v b ) « v d - v b ~ A. 


(18) 


Eq. [l 8 ]says that the average cell size is a constant roughly 
independent of the growth rate, a prediction contradict¬ 
ing the well-tested exponential scaling with growth rate 
for the model organisms mentioned above. 

Furthermore, the above reasoning assumes that the 
master accumulation model can maintain a balanced cell 
cycle. But the master accumulation model cannot ro¬ 
bustly maintain a balanced cell cycle in face of noise as 


Eq. 14 demonstrates. Specifically, if a cell initiated repli¬ 
cation at volume Vi = A, then the duration to the next 
initiation event is 


U(0) 


T ln f 1 + A 

A V A 


(19) 


But the cell may proceed to initiate an extra round of 
replication if the noise is negative enough, £/r < [(C + 
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D)/t\ — (C + D)/t , which corresponds to an extra round 
of initiation at volume v\ = 2A — 5v. The next initiation 
will then occur after 

< 20 > “ Ml + 2A^) S ‘» & (J) T ' (20) 

Implying that: 



-0.3 


0.0 0.5 1.0 1.5 2.0 

a C+Dl a T 



0.0 0.5 1.0 1.5 2.0 

a C+Dl a T 


U (20) < u (O) (21) 

Eq. [21] says that cells with more origins will initiate 
faster than those with less, giving rise to cells with av¬ 
erage durations between birth and division not equal to 
r. In other words, the master accumulation model does 
not robustly regulate the initiation process to maintain 
a balanced cell cycle. Indeed, simulations of the master 
accumulation model do not converge to a balanced cell 
cycle. Likewise, |5j carried out simulations of the master 
accumulation model and obtained ’’widely abnormal cell 
size distributions.” Given the above inconsistent predic¬ 
tions, the master accumulation model can be ruled out as 
a possible regulation strategy for replication initiation. 


FIG. 6. Pearson (left) and Kendall (right) correlations be¬ 
tween cell sizes at birth Vb and cell size differences between 
birth and division Av against (jc+d/^t suggests that varia¬ 
tions in C + D are small. 

in durations between initiation and division, but these 
variations are not accounted for by the multiple origins 
accumulation model. Variations of this nature give rise 
to the negative correlations between Vb and Ate 

Multiple origins accumulation predicts 
proportionality between cell size and the number of 
origins per cell 


Multiple origins accumulation suggests that 
variations in C + D are small 

In claiming that the master accumulation model gives 
incorrect correlations between growth rate dependent 
variables, [5] simulated the master accumulation model 
and reported negative correlations between cell size at 
birth Vb and cell size differences between birth and di¬ 
vision Av, whereas none is observed experimentally. In 
contrast to claims in [9], the negative correlations do not 
provide evidence against the multiple origins accumula¬ 
tion model nor the master accumulation model. Instead, 
the negative correlations provide evidence that variabil¬ 
ity in the durations from initiation to division C + D 
should be small. Indeed, the multiple origins accumula¬ 
tion model, because of its reduction to the incremental 
model of size control, predicts no correlations between 
Vb and Av, given that variations in C + D are small 
compared to variations in r. Simulations assuming that 
durations from initiation to the corresponding division 
are normally distributed with mean C + D and standard 
deviation ac+D show that the correlations between Vb 
and Av become increasingly negative as ac+o/a T in¬ 
creases, as seen in Figure [6] Figure [6] shows that as long 
as ac+ d j uy < 0.3, the correlations between Vb and Av 
will be close to zero. This is intuitive because when ac+D 
is small compared to ay, fluctuations in cell sizes at birth 
arise due to variations in cell sizes at initiation, but these 
variations are negatively fed back into the multiple ori¬ 
gins accumulation model as explained in Section [] On 
the other hand, when ac+D is comparable to ay, some 
fluctuations in cell sizes at birth arise due to variations 


The simultaneous regulation of cell size and the num¬ 
ber of origins per cell in the multiple origins accumulation 
model gives rise to a strict relationship between the two 
variables. In the multiple origins accumulation model, 
the average cell volume at birth Eq. 12 is exponentially 
dependent on the growth rate, while the average number 
of origins per cell Eq. [3] also scales exponentially with 
the growth rate. Therefore we have that 


( v b ) ~ A (O) 


( 22 ) 


That is, the multiple origins accumulation model pre¬ 
dicts that given a fixed volume increment per origin A, 
the average volume at birth and the average number of 
origins per cell will scale appropriately with respect to a 
varying [C + D )/r to give rise to the above approximate 
proportionality. The critical size regulation strategy pro¬ 
posed by Donachie assumed this proportionality 0 > but 
is inconsistent with measured correlations in variables in 
E. coli because thresholding size at any point in the cell 
cycle washes away the memory of the initial conditions, 
and therefore leads to a vanishing correlation coefficient 
between size at birth and size at division - contrary to 
measurements; on the other hand, multiple origins accu¬ 
mulation predicts this proportionality, and is consistent 
with measured correlations EH3- Simulations of the mul¬ 
tiple origins accumulation model confirm that cell size 
is indeed approximately proportional to the number of 
origins per cell, Figure [7] We emphasize that the ap¬ 
proximate proportionality is a property predicted by the 
multiple origins accumulation model. In contrast, other 
strategies that do not regulate the number of origins, such 
as the master accumulation model, would not predict it. 
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FIG. 7. (O) against (vt) shows proportionality between 

cell size and the number of origins per cell. Simulations are 
as in Figure [3j with a varying t = 20 to 100 mins, a fixed 
C + D = 70 mins, and o>/t = 0.2. Dashed line plots Eq. [3] 
and Eq. [12] 


Multiple origins accumulation predicts bimodal cell 
sizes at initiation 



FIG. 8. Distribution of volumes at initiation is bimodal. 
Simulations are as in Figure [3] with r = 100 mins, C + D = 
70 mins, and o> = 20 mins as a specific, slow growth example. 
Square symbols separate volumes at initiation into two sub¬ 
populations, one whose cells initiated the predicted number 
of rounds of replication at volumes near JoA, and another 
whose cells initiated extra rounds of replication near twice 
that volume 2/oA. Solid lines plot lognormal distributions 
with means 7 qA and 2/oA and variances a1. = 4a^/3r 2 [5]. 


In addition to the approximate proportionality be¬ 
tween cell size and the number of origins per cell, the 
multiple origins accumulation model predicts that the 
distribution of cell sizes at initiation will be approxi¬ 
mately bimodal because cells will initiate extra rounds 
of replication when the noise is large enough. Simula¬ 
tions show that the distribution of cell sizes at initiation 
is indeed bimodal, with one large peak corresponding to 
a subpopulation whose cells initiated the expected num¬ 
ber of rounds of replication and a smaller subpopulation 
whose cells initiated extra rounds of replication, as seen 
in Figure[8] Naively, because the distribution of cell sizes 
is lognormal in the multiple origins accumulation model, 
the distribution of cell sizes at initiation should be ap¬ 
proximately the sum of two lognormal distributions with 
means OqA and 20oA, where Oq = 2 ^ C+D ^ T ^, and the 
ratio between the frequencies of the two peaks equal to 
the probability that £/r < [(C + D)/t\ — (C + D)/t. 
However, the value OqA overestimates the average sizes 
at initiation of cells that initiated extra rounds of repli¬ 
cation because of correlations between the volumes at 
initiation and the probability for extra rounds of replica¬ 
tion. The correlations arise from Eq. [2] which says that a 
smaller volume at initiation correlates with a larger prob¬ 
ability for extra rounds of replication during the current 
cell cycle. The bimodal distribution of cell sizes at ini¬ 
tiation highlights how the multiple origins accumulation 


model, without invoking other mechanisms, can robustly 
maintain a balanced cell cycle despite noise in the initia¬ 
tion process and is an experimentally testable prediction 
of our model. 


Multiple origins accumulation predictions are 
consistent with experiments on mutants 

Experiments on mutants of E. coli and B. subtilis pro¬ 
duced results consistent with the predictions of the multi¬ 
ple origins accumulation model. Experiments have shown 
by manipulating the cell size of E. coli via mutations that 
a decrease in v is correlated with a decrease in C + D and 
that a decrease in v is also correlated with a decrease in 
(O) @nj. In the language of the multiple origins accumu¬ 
lation model, v is controlled by A and ( C + D)/t. We 
assume that the mutations left unchanged the replication 
initiation mechanism in E. coli , so that A is a constant 
throughout. But r remained approximately constant 
with changes in v, so a decrease in v must correspond to a 
decrease in C+D. One particular measurement reported 
for E. coli cells an average doubling time r ss 25 mins for 
wildtype and mutant cells, C + D ss 40 mins + 20 mins 
for wildtype cells, and (C + D) ~ 30 mins + 20 mins 
for mutant cells [20]. Given these values and Eq. 
the relative change in cell sizes corresponding to the re- 
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ported difference in the durations from initiation to di¬ 
vision should be 2^ C+D ' 1 ~ t - c+D '>)/ T « 0.75. That is, the 
mutant cells should be 0.75 times the size of the wild- 
type cells. This is in excellent agreement with the 25% 
decrease in volume reported. Moreover, the size at initia¬ 
tion did not change for these smaller mutants, consistent 
with our model which predicts the size at initiation to 
depend only on A and not on C + D. In contrast, in 
the case of B. subtilis, C + D in smaller mutant cells re¬ 
mained constant, suggesting that A is the quantity which 
was changed. Based on this interpretation, our model 
would predict that both size at initiation and at birth 
would change proportionally. Indeed, it was found that 
both the average mutant cell size and the mutant cell size 
at initiation both decreased by approximately 35% [20] . 
The above experiments also observed the predicted ap¬ 
proximate proportionality between cell size and the num¬ 
ber of origins per cell in both E. coli and B. subtilis [20] . 
It remains to be shown that cell sizes at initiation fall 
into an approximate bimodal distribution. The agree¬ 
ment between experimental results and the predictions 
made by the multiple origins accumulation model speaks 
to the importance of regulating simultaneously cell size 
and the number of origins. 

DISCUSSION 

The multiple origins accumulation model proposes that 
replication initiates upon the accumulation of a critical 
amount of initiators per origin. If the initiators are ex¬ 
pressed as in the autorepressor model, this strategy corre¬ 
sponds to Eq. |T] which in turn reduces to the incremental 
model of size control, which predicts distributions, cor¬ 
relations, and scalings consistent with existing measure¬ 
ments. Specifically, the average cell size scales exponen¬ 
tially with the growth rate of the cell, Eq. [12] as does 
the average number of origins per cell, Eq. [3] The model 
robustly regulates both cell size and the number of ori¬ 
gins per cell such that cell size is approximately propor¬ 
tional to the number of origins per cell, Eq. [22] These 
predictions are consistent with existing experiments on 
E. coli and B. subtilis [ 20] , A proportionality between 
ploidy and cell size has also been observed in other or¬ 
ganisms, including yeast J2Tj. The multiple origins accu¬ 
mulation model is a general regulation strategy that may 
illuminate the source of the approximate proportionality 
between cell size and the number of origins across organ¬ 
isms. 

An essential feature of the multiple origins accumu¬ 
lation model is the tight coupling between chromosome 
replication and cell division. The differences between the 
multiple origins accumulation model and the master ac¬ 
cumulation model emphasizes this coupling and the im¬ 
portance of regulating the timing of initiation. By nega¬ 
tively regulating cell size in response to the number of ori¬ 


gins via Eq. |TJ the multiple origins accumulation model is 
able to maintain a balanced cell cycle and achieve robust¬ 
ness in face of noise in the initiation process. However, 
the master accumulation model described by Eq. [13] is a 
regulation strategy without such a feedback mechanism. 
The master accumulation model is unable to maintain a 
balanced cell cycle and does not predict the exponential 
scaling of cell size. This suggests that regulation strate¬ 
gies must account for the number of origins per cell in 
order to regulate appropriately the frequency of division. 

The coupling between the number of origins to the di¬ 
vision frequency could be demonstrated via a shift-up 
experiment. It was found for E. coli that cells maintain 
their rate of division for a duration of C + D after a shift- 
up, a phenomenon known as rate maintenance [221 251 . 
The multiple origins accumulation model naturally ac¬ 
counts for rate maintenance, because division always oc¬ 
curs at time C + D after initiation (chromosome replica¬ 
tion rate is independent of growth rate). Furthermore, 
the model offers a robust mechanism to regulate, after 
a transient, the number of origins per cell appropriately 
with the new growth rate via Eq. [3] The existence of a 
rate maintenance period implies that division is coupled 
to replication initiation, and the incremental model ap¬ 
plied at birth and division is a valid phenomenological 
description only at stationarity. Instead, it is the under¬ 
lying molecular mechanism of replication initiation that 
dictates the frequency of division. 

Although the multiple origins accumulation model cap¬ 
tures many aspects of the coupling between replica¬ 
tion and division, experiments with minichromosomes 
suggest that the molecular mechanism is more compli¬ 
cated. Minichromosomes are plasmids containing the 
oriC sequence coding for chromosomal origins. In gen¬ 
eral, minichromosomes initiate replications in coordina¬ 
tion with chromosomes and do not affect the growth 
properties of the cell, such as the doubling time or the av¬ 
erage cell size [24]. However, if more than ~ 40 minichro¬ 
mosomes are present in a cell, replication initiation is no 
longer synchronous, the doubling time increases, the av¬ 
erage number of origins per cell decreases, and the aver¬ 
age cell size decreases [25] . Another experiment inserted 
a second origin into E. coli chromosome and observed 
again that the extra origin does not affect the growth 
properties of the cell (25] ■ These result points to a more 
complicated molecular mechanism than accumulation of 
initiators per origin. Several mechanisms have been sug¬ 
gested, but none has been completely satisfactory. For 
example, the master accumulation mechanism discussed 
in Section]] is ruled out for being unable to robustly regu¬ 
late the number of origins. Another plausible regulation 
strategy is one in which replication initiates when a criti¬ 
cal ratio of active to inactive initiators is reached [4]. The 
validity of this strategy remains to be tested. 

The molecular mechanism underlying the regulation of 
replication initiation is yet to be unraveled, but here we 
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have given significant constraints regarding the potential 
mechanisms. Specifically, this work and previous works 
have shown that the molecular mechanism in question 
should satisfy both the incremental model of size con¬ 
trol and the mathematical form of the multiple origins 
accumulation model described by Eq. [T] so that the pre¬ 
dicted distributions, correlations, and scalings remain in¬ 
tact and consistent with existing experiments. 
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APPENDIX: DERIVATIONS OF THE AVERAGE 
NUMBER OF ORIGINS PER CELL 

The average number of origins per cell has been calcu¬ 
lated previously in two distinct derivations [2 CH|. The 
two derivations seemed to make different assumptions to 
arrive at the same conclusions, bringing into question the 
necessity of the underlying assumptions. Here, we repro¬ 
duce the two derivations and show that both derivations 
in fact make the same assumptions. 

The model of the cell cycle under consideration is due 
to Cooper and Helmstetter p]. In this model, replication 
initiation occurs on average every doubling time r. An 
initiation event then triggers a division event after a con¬ 
stant duration C+D , where C and D are respectively the 
constant duration required to replicate the chromosome 
and the constant duration between replication termina¬ 
tion and division. Given C', D , and r, we want to find the 
average number of origins per cell. The average number 
of origins per cell is defined as ( O) = (OtotaiMO, where 
Ototai is the total number of origins in a population of 
cells, N is the number of cells in that population, and 
brackets denote the ensemble average. 

First, we reproduce the derivation due to Cooper and 
Helmstetter [1], To calculate the average number of ori¬ 
gins per cell, we must first define the probability dis¬ 
tribution underlying the ensemble average. In an asyn¬ 
chronous population of exponentially growing cells, the 
cells must be exponentially distributed in the cell cycle 
for ensemble averages to be stationary with respect to 
time [28]. Defining cell age a = 0 at birth and a = 1 at 
division, the exponential distribution of cell ages is 

p (a) = (ln2) 2 1 ~ a . (23) 

We can now calculate the desired ensemble averages. For 
example, if 0 < C + D < 1, then a cell younger than 
(r — (C + D)) /r will not be replicating its chromosome 
and will have only one origin, whereas a cell older than 
(r — (C + D)) /r will be replicating its chromosome and 
will have two origins. We have assumed that the amount 
of time a cell spends with more than two origins is neg¬ 
ligible, which is plausible for weak noise in the initiation 
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process. The average number of origins per cell in an 
asynchronous, exponentially growing population is then 


(O) = In 2 



2 - 2 (C-kd)/t ^ + 2 ^2 (c+d,/t - l) 


= 2 (C+d)/t 


(24) 


Similarly if 1 < C + D < 2, then a cell must initiate 
replication not only for its daughter cells, but also for 
its granddaughter cells. In this case, a cell with cell 
age less than (2r — (C + D)) It will not be replicating 
its chromosome for its granddaughters and will have 
only two origins, whereas a cell with cell age more than 
(2r — (G + D)) /r will be replicating its chromosome for 
its granddaughters and will have four origins. Again, we 
have assumed weak noise. 

Thus, we see that 


/■AT 

p (O = O 0 ) = ln2 / 2 1 ~ a da, (25) 

Jo 

p{0 = 2O 0 ) = In 2 [ 2 1 ~ a da, (26) 

J AT 

where O 0 = 2^ c+D ^ T i, and AT 
(L(C + D)/t\ + 1) t — (C + D). Simplification gives 


(27) 


which generalizes Eq. 24 and is valid for any C + D 
and r. The two assumptions made in this derivation are 
that the population is growing exponentially and that the 
population has reached a stationary distribution of cell 
ages. 

Next, we reproduce the derivation due to Bremer and 
Churchward E3- Assuming exponential growth, the 
number of cells must grow exponentially as N oc 2 4 / T . 
Similarly, the total number of origins must grow at the 
same exponential rate so that O tot oc 2 4 / T . But an ini¬ 
tiation event triggers a division event after a constant 
duration C + D, so the number of cells must on aver¬ 
age lag behind the total number of origins by 2^ C+D ^ T . 
The average number of origins per cell must then be 
(O) = ( Otot/N) = 2( c+d ' ) / t . Although the distribution 
of cell ages was not explicitly involved in this derivation, 
the assumption of a stationary ensemble average under 
exponential growth is satisfied if and only if the distribu¬ 
tion of cell ages is exponential (28]. 

The exponential distribution of cell ages is not always 
realized in experimental setups. For example, single-cell 
experiments that track a lineage of cells, such as those 
in [5] , will follow a different distribution, as discussed in 
[29] . Experiments that track a single cell will follow a 
uniformly distributed cell age. In that case, Eq. [27] is 
replaced by 


(O) = o 0 ( 1 + £±5 _ l c±d 

T T 


(28) 


Simulations tracking a population of cells with uniformly 
distributed cell ages confirm this result. The differ¬ 
ences between Eq. [27] and Eq. [28] do not significantly 
change the predictions of the multiple origins accumula¬ 
tion model. 


(O) = 2^ c+d ^ t , 












